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We discuss the ten dimensional black holes made of DO-branes in the regime where 
the effective coupling is large, and yet the 11D geometry is unimportant. We suggest that 
these black holes can be interpreted as excitations over the threshold bound state. Thus, 
the entropy formula for the former is used to predict a scaling region of the wave function 
of the latter. The horizon radius and the mass gap predicted in this picture agree with the 
formulas derived from the classical geometry. 
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The 10 dimensional black hole consisting of DO-branes in the regime where the clas- 
sical 10 dimensional geometry defined by the string metric is reliable has been hard to 
understand. This difficulty is caused partly by the large effective coupling constant in the 
DO-brane quantum mechanics, and partly by the fact that unlike other higher dimensional 
D-branes, DO-branes are highly dynamical due to the Heisenberg uncertainty relation. 

Our basic observation in this paper is that the black hole in the regime of interest can 
be regarded as a small excitation over the large N threshold bound state. Thus, it is hope- 
less to make a first principle calculation concerning the black hole without understanding 
the bound state first. What we shall do here is to postulate that there is a scaling regime 
where the one-particle wave function obeys a power law, matching the entropy formulas 
fixes the exponent. As a consistency check, we show that the correct formula for the hori- 
zon radius emerges naturally in this picture. We also estimate the mass gap in the black 
hole background. 

The near horizon geometry of a DO-brane black hole is given by 

ds 2 = -fdt 2 + f~ l du 2 + Vxu- s / 2 dn 2 , 

A\ 3/4 (!) 
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where 
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We have set a' = 1. It was argued by Polchinski that there is a duality between M theory 
defined on the background geometry of the threshold bound state and quantum mechanics 
of DO-branes in the large N limit || . For related discussions, see @ . 

This duality reduces to that between the IIA theory and DO-brane quantum mechanics 
when the curvature is smaller than the string scale, and the dilaton is sufficiently small. 
The conditions are just 

1 « v/AtV 372 « N 2 ' 7 . (3) 

The Hawking temperature is 

T = -^d ugoo (U = U ) = ^\~ 1/2 U* /2 . (4) 
In the M theory unit, we have 



r = Uoa' = N 1 /%/\[3R)- 2/5 , 
1 
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where we ignored a numerical factor. Define the effective coupling by A e // = NR [3 /I, 
the entropy of the black hole is given by 
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S = N 2 X 




(6) 



again a numerical factor is dropped. Conditions (|3]) are rewritten in terms of the effective 
coupling 



Since the effective coupling is large, it is hopeless to describe this black hole using the 
perturbative quantum mechanics of DO-branes. Also notice that there is a upper bound 
on the effective coupling, this is to avoid the M theory region where the effective 11th 
dimension becomes larger than the Planck length. Thus, the 10 dimensional black hole 
we are aiming to describe in this paper is different from the matrix black holes discussed 
in This upper bound also indicates that for a given temperature, we do not have to 

hold the 't Hooft coupling fixed. Indeed, the natural limit in matrix theory is to hold the 
Yang-Mills coupling fixed. It is unclear whether a different phase should occur when one 
crosses the upper bound for X e ff- The entropy of the black hole is of order iV 8 / 7 at this 
bound, this is still larger than iV where a localization transition occurs, for a systemacally 
discussion on the phase diagram, see ||. We shall here adopt a conservative viewpoint, 
and obey the upper bound. 

As in the background of any black hole in any dimension, there is a mass gap. The mass 
gap is determined by the Hawking temperature. A WKB solution of the dilaton equation of 
motion clearly indicates that indeed this mass gap is order of the temperature [0 . We shall 
not quote the precise formula here, since we will not try to calculate the exact coefficient. 
Again, this behavior is hard to understand in terms of the perturbative nonabelian quantum 
mechanics, which is valid only in the high temperature regime. A rather straightforward 
calculation shows that the one-loop mass gap m? is actually negative, and is proportional 
to A e //T 2 , see the appendix. 

Recently it has been argued that the size of the large N threshold bound state grows 
with N as N 1 / 3 @,|]. The bound state is a highly nonabelian object. It is impossible 
to separate the DO-brane degrees of freedom from the off-diagonal elements. It is still 
reasonable to talk about the density of DO-branes. The situation is similar to the matrix 
membrane 0, where the object is also highly nonabelian, nevertheless so long if N is 



1 < Xeff < N 
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sufficiently large, it is meaningful to define DO-brane or the longitudinal momentum density 
along the membrane. 

We shall not repeat the argument which leads to the estimate of the size of the bound 
state L ~ N x ^l p . For our purpose, it is sufficient to know the order of magnitude of a few 
relevant quantities. The size can be measured by the quantity 1 / N (trX 2 ) = L 2 , X is one 
of the nine bosonic matrices. The typical frequency of an off-diagonal element is of order 

i 

u = N^Rl' 2 , (8) 

thus the "zero-point energy" of an off-diagonal state is of this order. This huge amount 
of energy is supposed to be canceled by contribution from the fermionic fluctuations. If 
we assume that the black hole is a "small excitation" over the threshold bound state, the 
surplus energy must come from excitations of the off-diagonal elements. This is because 
the energy is of order iV 2 when \ e jf is held fixed. The typical energy carried by such 
elementary quanta is just T. This is much smaller than u since u/T = A^jj ^> 1. Indeed 
this is a small fluctuation over the bound state configuration as long as the bosonic degrees 
are concerned. So the wave function will not be significantly changed by the fluctuation. 
This is our basic observation. 

It is easy to see that the condition X e ff ^> 1 is also equivalent to the condition r$ <C L. 
Physically, one would imagine that it is sufficient to know the shape of the wave function 
in this region in order to calculate the thermodynamic quantities of the black hole. This 
will turn out to be true. For our purpose, it it enough to know the density of DO-branes 
in this region, p(r) ~ ^(r)! 2 . We shall adopt the following scaling ansatz 

p(r) = CNL- Q ^)\ (9) 

where a is to be determined, C is a numeric constant which we will drop henceforth. We 
used rotational invariance in this ansatz. Obviously, L must be the only relevant scale in 
this problem. Beyond L, there may be another scaling region which can not be probed by 
black hole. This is the perturbative region in quantum mechanics, however. 

The above ansatz is valid only when we smear over scales much larger than the average 
spacing between two neighboring DO-branes. As we have emphasized, it is difficult to define 
such spacing with a highly nonabelian configuration. When we make an observation only 
in one direction, say X, it makes sense to define such a spacing. It is L/N ~ N~ 2 ^ s l p: a 
rather small quantity in the large N limit. The high density in one dimension is caused by 
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the projection. A naive definition of the true spacing is a ~ (L 9 /iV) 1 / 9 ~ N 2 / 9 l p . This is 
a rather large quantity. Lacking a better definition, we shall use this naive one here. 

The frequency for a diagonal element stretched between a DO-brane at point X and 
another DO-brane at point Y is i?|X — Y|Z~ 3 , thus the Boltzmann weight is a function 

/(/3i?|X-Y|/- 3 |). (10) 

We shall not specify this function. For a bosonic degree of freedom, f(x) = (e x — 1) , and 
for a fermionic degree of freedom, it is f(x) = (e x + l) -1 . To include all possible states, 
/ can be more complicated. (There is no chemical potential, since individual elementary 
quanta associated to an off-diagonal element are not conserved.) The contribution of these 
states to the free energy is another function g(x). For a boson, it is ln(l — e~ x ) and for a 
fermion, it is ln(l + e~ x ). The total free energy is 

(3F = J p(X)p(Y)d 9 Xd 9 Yg(f3R\X - Y|/ p " 3 |). (11) 

Rescaling X — > LX and using the ansatz ([|) we obtain 

(3F = N 2 J d 9 Xd 9 Y\X.\ a \Y\ a g(\ 1 e f f \X - Y|). (12) 

The entropy can be obtained using the standard formula S = (3E — (3F. It is obvious 
that the entropy is a function of N 2 and A e // only, just like j3F. This result justifies our 
observation that indeed the black hole is a thermal excitation over the threshold bound 
state. 

— 1/3 

To obtain the right dependence on A e //, let us do a further rescaling X — > A e j^ X. 
The free energy reads 

0F = N 2 \; f (2 f a+18)/3 J d 9 Xd 9 Y\X\ a \Y\ a g(\y. - Y|). (13) 

Now the integral is a pure number. Matching onto the scaling in (^|) determines the 
exponent 

« = --. (14) 

10 V ! 

This is our main result in this paper. With this exponent, it is obvious that the integral 

J d 9 X p(X) is convergent at the origin, and divergent at large X. This is just fine, since 

our ansatz (|9|) is supposed to be valid only for |X| < L. We shall see shortly that this 

restriction agrees with the black hole physics as governed by formula (|i~3l) . Beyond the size 



L, we enter into the perturbative region of quantum mechanics. For a threshold bound 
state, the decay of the wave function at sufficiently large distance also obeys a power law. 
The exponent in the asymptotic region must be smaller than —9. 

We still have to worry about the convergence of the integral in (|i~3|) . The function g(x) 
damps fast for large x, peaks at x = 0. Integrating over Y first, we shall obtain another 
factor |X| a for large |X|, thus integration over X is convergent for large |X|. For small X, 
integrating over Y first will give rise to a order 0(1) number, so integration over X near 
the origin is also convergent. We conclude that the integral in (^) is a finite pure number. 

i . 1/3 

The contribution to the integral in ( p~3|) mainly comes from the region |X| ~ ^^ e ff > 
in terms of the coordinates before rescaling. This size is smaller than L, but independent 
of N. Thus it can not be identified with the horizon radius (|j). There may be a few 
different definitions of the horizon radius. One simple definition is the subtracted quantity 
l/iV(trX 2 ). By subtraction we mean that the physical contribution to the average comes 
only from the thermal fluctuations. This is reasonable, since by scattering an object against 
a threshold bound state, the wave function size should be not probed, as a basic assumption 



in a holographic theory. Another more precise definition is proposed in [|10| . The proposal 
is that the horizon radius is where a tachyonic mode develops between a probing DO-brane 
and the black hole. The authors of |I(J argue that the horizon radius is determined by the 
first negative eigen-value of the following matrix 

(M 2 j ) = (5 ij X 2 + 2[X i ,X j }). (15) 

It is expected that the negative eigen-value will be of order l/iV(trX 2 ). Since we are not 
trying to compute the exact coefficient, we will not endeavor to use the above definition. 

The amplitude X 2 of an off-diagonal element stretched between point X and point Y 
is nlpl'K — Y| _1 , when it is in its n-th harmonic state. The thermal average is than 

-±-ftfR\X-Y\l-*), (16) 

where / is the Boltzmann factor. We have the following thermal average 

1_, I 3 
N 

Performing the same rescaling as before, we obtain 
i A/7 3 f 

i-{trX 2 )f3 = -^Xjl 15 / d 9 Xd 9 F|X|- 81 / 10 |Y|- 81 / 10 |X - Y|" V(|X - Y|). (18) 
iV L 31 J 



■(trX 2 ) p = i J rf 9 rfy P (X)p(Y)|X-Y|-V(/3^|X-Y|/ p - 3 ). (17) 



Now, the integral is also convergent, since the superficial pole in |X — Y| _1 as well as in 
/(|X — Y|) is not severe in 9 dimensions. Thus this integral is a pure number. Substituting 
what we know about L and X e ff, the above quantity is of order 

N 2 'Hf'\f3R)-^\ (19) 

This is precisely r%, as in (|5|). ro can be also expressed as LA e // 15 . There are three scales 
appearing in our problem 

l > la;// 15 > la;// 3 . (20) 

We are fortunate that it is not the last scale being identified with the horizon radius. The 
correct horizon radius we obtained is a strong consistency check of our picture. 

The next thing we wish to do is to identify the mass gap. Our definition of the mass 
gap is the quantity m 2 appearing in the following effective action 

S = I J dtix ((A) 2 + m 2 X 2 + ...), (21) 

where we use the Euclidean time, the dots denote high order terms. The quadratic term 
trA 2 comes from thermal fluctuations. We expect the answer be m 2 ~ T 2 . As explained 
in the appendix, this is going to be a nonperturbative result. 

It is rather easy to see how to compute m 2 in our picture. It is just the average mass 
square of stretched strings. From the quadratic term Rl~ 6 tr[Xi, Xj] in the matrix action, 
we find the right quantity to be averaged over is R 2 l~ 6 \K - Y| 2 . An additional factor R 
comes from our definition in (pi|). Obviously, up to a numerical factor, one can replace 

1 /3 

|X - Y| 2 by (L\~ f y) 2 . Therefore the mass gap is 

m 2 ~ J R 2 /; 6 (LA;// 3 ) 2 = T 2 . (22) 

We see that all factors depending on R, l p and N cancel. In a way, we should not be 
surprised by this result. 



It is possible to follow the strategy of |Tl] to compute the Hawking radiation rate. It 



may also be possible to use the exact coefficient in the entropy formula to fix the coefficient 
in the formula of DO-brane density, then coefficients in other physical quantities may be 
computed within our picture. We are only beginning to investigate the fascinating subject 
of the large N DO-brane threshold bound state, and its role in the duality between the large 
N DO-brane quantum mechanics and M theory. We have only performed some qualitative 
analysis. Hard work remains to be done. 
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I understand that D. Kabat and G. Lifschytz have been investigating the ten dimen- 
sional black holes along a different line. 

Acknowledgments This work was supported by DOE grant DE-FG02-90ER-40560 and 
NSF grant PHY 91-23780. Part of this work was carried out during a visit at Rutgers 
University and a visit at KIAS at Seoul. Both institutions are gratefully acknowledged for 
their warm hospitality. I have benefited during the course of this work from conversations 
with T. Banks, M. Douglas and P. Yi. 
Appendix 

There are a few papers discussing DO-brane statistical mechanics in the perturbative 



regime [12]. To our knowledge, the one-loop mass gap has never been computed. We will 



give a formula at the one-loop level here. 

The Euclidean action of DO-brane quantum mechanics is 

Se = [ (jR iD ^ - ^ X " X ^ + 9D ^ + ^ X ^) ' (23) 
Performing rescaling r — > /3r, Xi — > \fj3RXi, we reach the following action 

S E = J drtr Q(^rX) 2 - £[Xi, X 3 } 2 + 6D T 6 + gOy^, 6]\ , (24) 

where the effective coupling is g 2 = ((3R) s I~ 6 (4it 2 )~ 1 . 

The mass gap is the coefficient appearing in the effective quadratic potential. One 
way to compute the static potential is to use the background field method. The first step 
is to fix the gauge. In general, the gauge field A can be put into a diagonal form, each 
diagonal element represents a U(l) holonomy. We are only interested in the mass gap, so 
we set A = 0. For fermions, since the action is quadratic, one can even compute the exact 
static potential coming from this part. The simplest way to compute it is to start with 
the Hamiltonian formalism. For a given static background Xi, the Hamiltonian is simply 

H, = |tr {Oyi[X it 6]) = \e a a N aa , b pd% (25) 

where a, b are indices of the su(N) Lie algebra, and a, (3 are spinor indices. We have rescaled 
0's such that {6%, 6 b p } = 25 ab 5 a p. The 16(iV 2 - 1) x 16(iV 2 - 1) matrix J\f is 

K a , b ? = \igifxtf ahc . (26) 



It is easy to check that Af is an antisymmetric matrix. 

Now Q a a form a 16(iV 2 — 1) dimensional Clifford algebra. The fermionic Hilbert space 
is a spinor representation of this algebra. The static potential is encoded in tre - ^. Since 
Af is an even rank antisymmetric matrix, by an orthogonal rotation, it can be put into the 
Jordan form. Once this done, we see that 

tre"^ = (det(e^ + e"^)) V " . (27) 

The determinant is invariant under the orthogonal rotation, thus the above formula is valid 
for a general N '. For a small Af, one can perform a perturbative expansion. The first term 
in the effective action is quadratic in Af. Indeed we have 

where the trace is taken over (a, a). In the end, we have 

tre""' =e ^ 2 trX 2 +..._ (28) 

Apparently, this expansion is valid when the coupling A e // = Ng 2 is small. We see that 
the contribution of fermionic fluctuations to the mass gap m 2 is negative. 

It is a little more complicated to compute the contributions from bosonic fluctuations. 
We separate bosonic matrices into a sum of the background and the fluctuation Xi + Y t . 
The one-loop contribution is determined by the quadratic part in the fluctuation Yi in the 
action 

Sb = \{d T Y?) 2 + ±Y?Mai,bjY}, (29) 

with 

M a i,bj = g 2 [-i[Xi, Xj] c f abc + Xj:X^f ace f hde 6ij - XjXff ace f bde ]- (30) 

The path integral can be performed in two steps. First one converts it into the form 
tre~ Hb . This gives the result 

tre~ Hb = (det[sinh(VA4/2)]) 1 . (31) 

This result is singular when A4 = 0. The problem is caused by zero modes. We have 
already included the zero modes in the background fields, the contribution of this part 
must be subtracted. The subtracted result is then 



For small Ai, the contribution to the effective action from bosonic fluctuations is then 



-L t rM = -^NghrX 2 + .... (33) 

Although this contribution to the mass square is positive, its magnitude is smaller than 
that from fermionic fluctuations. The overall mass gap m 2 is negative. If we rescale back 
to the original time and fields X, we see that m 2 ~ Ng 2 T 2 ~ A e //T 2 . 

It is not surprising to have a negative mass square at a finite temperature. The pres- 
ence of temperature breaks supersymmetry, thus the cancellation between bosonic fluctu- 
ations and fermionic fluctuations is spoiled. The phase space of bosons is overwhelmed by 
that of fermions, that at short distance there is a repulsive force. This agrees with the 



result for two DO-branes in the last paper of |12|], where integration over the holonomy of 
the gauge field is also taken into account. We expect that high order terms in the static 
potential takes over for large X, thus the net force will be attractive. 

As we have seen in the main body of this paper, the mass gap m 2 is positive in the 
strong coupling regime, and is independent of X e ff- When one tunes down X e ff, there 
must be a point where the mass gap vanishes. One might take this point as a phase 



transition point. This phase transition is generally predicted in 13 
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